Free Energy Landscape Of Simple Liquids Near The 
Glass Transition 



Chandan Dasgupta|[| and Oriol T. VallsJ 

f Centre for Condensed Matter Theory, Department of Physics, Indian Institute 
of Science, Bangalore 560012, India 

| School of Physics and Astronomy and Minnesota Supercomputer Institute, 
University of Minnesota, Minneapolis, Minnesota 55455, USA 

Abstract. Properties of the free energy landscape in phase space of a dense 
hard sphere system characterized by a discrctizcd free energy functional of the 
Ramakrishnan-Yussouff form are investigated numerically. A considerable number 
of glassy local minima of the free energy are located and the distribution of an 
appropriately defined "overlap" between minima is calculated. The process of 
transition from the basin of attraction of a minimum to that of another one 
is studied using a new "microcanonical" Monte Carlo procedure, leading to a 
determination of the effective height of free energy barriers that separate different 
glassy minima. The general appearance of the free energy landscape resembles 
that of a putting green: deep minima separated by a fairly flat structure. The 
growth of the effective free-energy barriers with increasing density is consistent 
with the Vogel-Fulcher law, and this growth is primarily driven by an entropic 
mechanism. 
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1. Introduction 

A liquid quickly cooled to temperatures below its freezing point enters a metastable 
supercooled state. At lower temperatures, the supercooled liquid undergoes a glass 
transition to a state in which it resembles a disordered solid. The dynamics of 
supercooled liquids near the glass transition exhibits ||] multi-stage, non-exponential 
decay of fluctuations and a rapid growth of relaxation times, features which are not 
fully understood. 

An intuitive description that is often used || [| for a qualitative understanding of 
the observed behavior near the glass transition is based on the "free energy landscape" 
paradigm. This description starts from a functional that expresses the free energy of 
a liquid in terms of the time-averaged local number density. At high temperatures 
(or at low densities in systems, such as hard spheres, where the density is the control 
parameter), this functional is believed to have only one minimum, that representing 
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the uniform liquid state. As the temperature is decreased to near the crystallization 
point, a new minimum representing the crystal, with a periodic modulation of the 
local density, should also develop. In the "free energy landscape" picture, a large 
number of "glassy" local minima of the free energy, characterized by inhomogeneous, 
aperiodic density distributions, also appear at temperatures below the equilibrium 
freezing point. If the system gets trapped in one of these glassy local minima as 
it is cooled rapidly, crystallization can not occur and the subsequent dynamics is 
governed by thermally activated transitions among some of the many metastable glassy 
minima. If the system visits many of these minima during its evolution over a certain 
observation time, it behaves like a liquid over such time scales: the time-averaged local 
density remains uniform. However, the dynamics in this regime, governed by thermally 
activated transitions, is slow and complex. In this picture, the glass transition occurs 
when the time scale of transitions among the glassy minima becomes so long that the 
system is confined in a single "valley" of the landscape over experimentally accessible 
time scales. The features of such a free energy landscape would be very similar to those 
found j^, H [7) in certain generalized spin glass models with infinite-range interactions, 
and in spin models ||, [| with complicated infinite-range interactions, but no quenched 
disorder. The behavior of these mean-field models exhibits a remarkable similarity 
with the phenomenology of the glass transition. These results suggest that the free 
energy landscape paradigm indeed provides a good framework for the understanding 
of the properties of supercooled liquids near the glass transition. Such a description 
requires numerically obtained information about the topography of the free energy 
landscape of these liquids. 

We have carried out several numerical studies of a dense hard-sphere system, using 
a model free energy functional proposed by Ramakrishnan and Yussouff (RY) [[To] , a 
discretized version of which exhibits JIT| , fl2| a large number of glassy local minima 
at densities higher than the value at which equilibrium crystallization occurs. [The 
control parameter for a hard-sphere system is the dimensionless density n* = po<7 3 , 
where po is the average number density in the fluid phase and a is the hard-sphere 
diameter; increasing (decreasing) n* has the same effect as decreasing (increasing) 
the temperature where the temperature is the control parameter.] From numerical 
studies [[l3| [l4|, |l5| of the Langevin equations for this system, we found that the 
dynamics changes qualitatively at a "crossover" density near n* = 0.95. The dynamics 
of a system initially in the uniform liquid state remains governed by small fluctuations 
near the liquid free energy minimum when the density is lower than n*. For higher 
values of n* , the dynamics is governed by transitions among the glassy minima. The 
time scales for such transitions were estimated from a Monte Carlo (MC) method in 
Ref. [ fl6| and found to increase rapidly with density. 

Here we report results of additional numerical studies in which a new approach 
to the free energy landscape is used. We have developed and used a new MC 
procedure that enables us to study transitions between different glassy minima and 
thus investigate the topography of the free energy surface in phase space. We have 
located a large number of glassy minima of the free energy so as to yield their statistical 
properties as a function of density. The total number of glassy minima is found to 
remain nearly constant as the density is varied in the range 0.94 < n* < 1.06. The 
free energies of the glassy minima are distributed over a wide range between the 
free energy of the uniform liquid and that of the crystal. An appropriately defined 
"overlap" between different glassy minima is also found to exhibit a broad distribution. 
We have found pairs of glassy minima that differ from each other in the rearrangement 
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of a very small number of particles. The height of the free energy barrier that separates 
two such minima is quite small. Such pairs may be identified as "two-level systems" 
which are believed ]l7]] to exist in all glassy systems. 

Our computation of the probability of transition from a glassy minimum to the 
others as a function of the free energy increment (see below) and the MC "time" t 
leads us to define an effective barrier height that depends weakly on t. The growth of 
this effective barrier height with density is consistent with a Vogel-Fulcher form |Q 
for a hard-sphere system [|l9| . The dependence of the effective barrier height on t and 
the density indicates that the growth of the barrier height (and the consequent growth 
of the relaxation time) is primarily due to entropic effects arising from an increase in 
the difficulty of finding low free-energy paths (saddle points) that connect one glassy 
local minimum with the others. 



2. Model and Methods 



We characterize our system by a free energy functional F[p] of the form 1 1 



F[p] = Fi(po) + k B T 



dr{p(r) ln(p(r)/p ) - Sp{r)} 



- (1/2) jdvj dv'C{\v-v'\)5p{v)5 P {v') , (1) 

where Fi(po) is the free energy of the uniform liquid at density po, and Sp(r) = p(r)—po 
is the deviation of the density p at point r from po- We set Fi(po) = 0. In Eq. (jl|) , T is 
the temperature and C (r) the direct pair correlation function ppfl of the uniform liquid 
at density po, which we express in terms of n* = poo~ 3 (a is the hard-sphere diameter) 
by making use of the Percus-Yevick p0[ approximation. The direct pair correlation 
function p0| of simple model liquids characterized by an isotropic, short-range pair- 
potential with a strongly repulsive core (such as the Lennard-Jones potential) is very 
similar to that of the hard-sphere system at high densities. Therefore, we expect our 
results to apply, at least qualitatively, to such dense liquids. 

We discrctizc our system by introducing a cubic lattice of size L 3 and mesh 
constant h in which the variables pt, i — 1,L 3 , are defined as pi = p(ri)h 3 , where 
p(vi) is the density at mesh point i. The dimensionlcss free energy per particle f[p] is 
f[pj = (3F[p\/N where N — po(Lh) 3 — n*L 3 a 3 is the total number of particles in the 
simulation box, (3 = l/lksT) and a is the ratio h/a. 

Ideally, one would like to start the system in a known glassy local minimum of 
the free energy, and investigate the topography of the free energy surface near the 
starting point by allowing the system to evolve, and finding out which configurations 
it subsequently visits and where it ends up. A conventional Metropolis algorithm 
MC procedure Jl6| is inefficient at doing this because at the relatively high densities 
studied here, it would take a very long time for the system to move out of the basin of 
attraction of the initial minimum. To obviate this difficulty, we have devised what we 
call a "microcanonical" MC method. The algorithm is as follows: we choose a trial 
value of what we call the free energy increment, AF, or Af if we are dealing with the 
dimcnsionless version of F. Then, starting with initial conditions which correspond to 
a local free energy minimum, we sweep the sites i of the lattice sequentially. At each 
step and site, we pick another site j at random from the ones within a distance a from 
i. We then attempt to change the values of pt and pj top(pt+Pj) and (1 — p){pi +Pj), 
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where p is a random number distributed uniformly in [0, 1]. The attempted change is 
accepted, and this is the crucial point, if and only if the free energy after the change 
is less than F max = F + AF where F is the value of the free energy at the minimum 
where we start the computation. The simulation proceeds up to a maximum "time" , 
t m , measured in MC steps per site (MCS). We perform a sweep over a range of values 
of AF, with the same initial conditions. If AF is smaller than the height of the lowest 
free energy barrier between the starting minimum and any other "nearby" minima, 
the system will remain in the basin of attraction of the starting minimum. As we 
increase AF, there will eventually be one or more minima that the system can find 
within a "time" t < t m . These minima are separated from the initial minimum by free 
energy barriers of height less than AF. As AF is further increased, additional minima 
will be made accessible, and since additional paths will become available between the 
initial minimum and the minima already accessible at smaller values of AF, these 
minima may be reached in fewer MC steps. Clearly, if one obtains the information of 
near which minimum the system is, and how long it takes to get there, one can begin 
to map out the free energy landscape. 

To find out which basin of attraction the system is in at time t, we save the values 
of the variables pi at relatively frequent time intervals At. These configurations are 
then used as the inputs in a minimization procedure fll|| that determines which basin 
of attraction the system is in. The entire procedure is repeated a number of times 
(the "number of runs") and averaged over. We have carried out this procedure at 
densities in the range 0.94 < n* < 1.06. We did not consider densities lower than 0.94 



because previous studies 14, |l5| show that the dynamics of the system is governed 
by transitions among glassy local minima only at higher densities. Since the Percus- 
Yevick approximation becomes less accurate at relatively high densities p0[ , values of 
n* > 1.06 were not considered. 

We used two different sets of the sample size L and the mesh size h, in one case 
commensurate with a close-packed lattice, and in the other incommensurate. The 
computationally more intensive part of our simulations was carried out for systems 
of size L — 15 with periodic boundary conditions and mesh size h = ct/4.6. No 
crystalline minimum was found for these incommensurate values. The other portion 
of the computations was performed for systems with L = 12 and h = 0.25c These 
values are commensurate with a fee structure and a crystalline minimum is found 
at sufficiently high densities. Because of the smaller size of these samples, we were 
able to explore more extensively several aspects of the problem under consideration. 
Our computations for the L = 15 sample were carried out for t m = 15000 MCS and 
At = 5000 MCS, whereas computations for the system with L = 12 were carried out 
to t m = 8000 MCS with At = 2000 MCS. A detailed discussion of how the glassy 
minima used in our study were chosen and the structure of these glassy minima may 
be found in Refs. @ H ||. 



3. Results 



During the evolution of the system, we monitor /3F and the maximum and minimum 
values of the variables pi,i = 1,L 3 . If the system fluctuates near one of the 
inhomogeneous minima, then the maximum value of pi would be much higher than the 
value (close to poh 3 ) it would have in the vicinity of the uniform liquid minimum. The 
system does not move to the neighborhood of the liquid minimum for the values of 
AF considered here. The total free energy remains nearly constant at a value slightly 
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lower than the maximum allowed value, F max = F + AF. 

In our analysis of the process of transition of the system from the initial 
glassy minimum to the basins of other minima, we define a "critical" value, A/ C (i) 
(or AF c (t)), of the free-energy increment Af (or AF) as follows: at every time 
investigated (i.e. times 5000, 10000 and 15000 MCS for L — 15 and times 2000, 
4000, 6000 and 8000 MCS for the L = 12 samples), we test, for increasing values 
of Af, what is the probability, P(Af,t), that the system has moved to the basin of 
attraction of a free energy minimum distinct from the starting one. This probability, 
which we obtain by averaging over a sufficient number (ten to fifteen) of runs, is (at 
constant time) zero for very small Af and rises toward unity as Af increases. At 
a constant Af, it increases somewhat with MC time, as the system explores further 
regions of phase space. We define Af c (t) as the value of Af for which, at that time, 
the switching probability reaches 1/2. Of course, P and A/ c are also functions of 
n*. The minima to which the system moves for values of Af close to or higher than 
Af c are, in general, different for different runs. This suggests that A/ c represents a 
measure of the free energy increment for which a relatively large region of phase space 
becomes accessible to the system. The system almost never returns to the basin of 
attraction of the initial minimum: after having left the initial minimum, the system 
cannot find its way back. 




Figure 1. (a) Example of the determination of the "critical" value A/ c , defined 
as the value of the free energy increment Af at which the transition probability 
P is 1/2. The black dots mark the intersections of the plots with the line P = 0.5 
(see text for a complete discussion). The data shown are for a sample of size 
L = 15, and three values (5000, 10000 and 15000 MCS) of the Monte Carlo time 
t. (b) Results for A/ c as a function of t at four different densities for a L = 15 
minimum. Results for other times and densities studied interpolate smoothly with 
the results shown. 

The procedure for determining A/ c is illustrated in Fig. [l]a, where we have shown 
the results for the transition probability P as a function of the free energy increment 
Af for a L — 15 minimum at n* = 0.99. It is clear from our data that the uncertainty 
in the estimated values of A/ c is ~ 0.05, the spacing between successive values of Af 
in the simulation. Typical results for Af c are shown in Fig. |l|b for the same L = 15 
minimum and four values of n*. Clearly, Af c is a weak function of t, and a stronger 
function of n*. The dependence of A/ c on t for fixed n* becomes more pronounced as 
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n* is increased. These dependences are analyzed later in this section. 

The large number of minimization runs we carried out locate a large fraction of 
the full collection of glassy minima of the free energy. For the "incommensurate" 
L = 15 sample used in our work, the number of minima we have located at each 
density is in the range of four to six. The "commensurate" L — 12 sample exhibits 
a substantially larger number of minima, one of which is crystalline (fee). A similar 
sensitivity of the number of local minima to the sample size and boundary conditions 
has been found in numerical studies |2^, of the potential energy landscape of model 
liquids described by simple Hamiltonians. In the following discussion of the statistical 
properties of the collection of glassy minima, we consider chiefly the results obtained 
for L — 12, for which we can produce significant statistics. 

The number of glassy local minima of the L — 12 system remains nearly constant 
as the density is varied in the range 0.96 < n* < 1.06. This number is close to 25. 
There is no systematic trend in the dependence of this number on the density. The 
free energies of these minima are distributed in a band that lies between the free 
energy of the uniform liquid (zero) and that of the crystal. The width of this band 
increases with n*. Since the number of minima is approximately independent of the 
density, this implies that the "density of states" of the glassy minima decreases as 
n* is increased. Let p(f3F)5 be the probability of finding a glassy minimum with 
dimcnsionless free energy between (3F — 6/2 and (3F + 5/2. We have calculated this 
quantity at different values of n* . Representative results at two densities, n* — 0.96 
and n* = 1.02, are shown in Fig. ||a. The values of 5 used are 4.0 and 8.0 for n* = 0.96 




Figure 2. (a) The "density of states" for glassy free energy minima, defined as 
the probability of finding a glassy minimum with free energy in a given range (see 
text). Results for L = 12 samples at two densities are shown, (b) The distribution 
P(q) of the overlap q defined in Eq. (j^J) . The distributions of the self-overlap and 
the overlap between different minima are shown separately for L = 12 samples at 
density n* = 1.02. 

and n* = 1.02, respectively. The range of (3F over which p((3F) is nonzero is clearly 
wider at the higher density. The consequent decrease in the values of p((3F) with 
increasing density is also clearly seen. Both distributions show peaks near the upper 
end, and tails extending to substantially lower values. However, the lowest free energy 
of the glassy minima is substantially higher than the free energy of the crystalline 
minimum. If the probability of finding the system in a glassy minimum is assumed 
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to be proportional to the Boltzmann factor , then only those minima with free 
energies lying near the lower end of the band would be relevant in determining the 
equilibrium and dynamic properties of the system. Our results indicate that the 
number of such "relevant" minima decreases with increasing n* . We find correlations 
between the free energy of a glassy minimum and its structure, similar to those found 
in Ref. jl^j]. Minima with lower free energies have more "structure" (as indicated by 
e.g. the heights of the first and second peaks of the two-point correlation function of 
the local density) and higher average density than those with higher free energies. 

We have also studied how the distributions of the local density variables in two 
distinct glassy minima differ from one another. The degree of similarity between two 
minima may be quantified in terms of their "overlap" [[7| . For the discretized system 
considered here, the overlap 9(1,2) between two minima labeled "1" and "2" may be 
defined in the following way: 

q(l, 2) = - J— max{R} ]T[pf } - p m ][p { X ~ Pav]- (2) 

Pav . 

Here, p\ and p^ are the discretized densities at the two minima and p av is the 
average value of the pi, which is assumed to be the same for the two minima. R 
represents one of the 48 symmetry operations of the cubic mesh, plus all translations 
taking into account periodic boundary conditions, R(«) is the mesh point to which 
mesh point i is transformed under R, and max{R} means that the R that maximizes 
the quantity on the right is to be taken. In Fig. |^b, we display the results for the 
distribution P(q) of q at n* = 1.02. The distribution of the self-overlap, defined as 
q(i,i) for the ith minimum, is also shown. As expected, the distribution of the self- 
overlap exhibits a sharp peak at a large value of q. The overlap between different 
minima exhibits a broad distribution with a peak at a small value of q, indicating that 
most of the glassy minima are rather different from one another. This distribution, 
however, extends to values of q as large as 0.6, indicating that there are a few pairs 
of glassy minima which are very similar to each other. For each value of n* , we 
find a small number (3-5) of such pairs of minima. The main difference between 
their structures comes from small displacements of just 2-3 particles. These pairs of 
minima are examples of "two-level systems" whose existence in glassy materials was 
postulated jl7j many years ago. 

We have also looked at how the quantity F c = Fq + AF C , varies from one minimum 
to another. While the free energy Fq of a glassy minimum varies over a wide range 
(see Fig. ^]a), the value of F c is nearly constant for each value of n* . This suggests a 
"putting green like" free energy landscape in which the local minima are like "holes" 
of varying depth in a nearly flat background. This structure also implies that there is 
a strong correlation between the depth of a minimum and the height of the barriers 
that separate it from the other minima: the barriers are higher for deeper minima. 

We now discuss the dependence of A/ c on the density n* and MC time t. Since 
the transition probability P is an increasing function of both A/ and t, Af c (n* ,t) 
decreases as t is increased (see Fig. [j]). In agreement with the previously observed [Q 
growth of the barrier-crossing time scale with n*, we find that A/ c is an increasing 
function of n*. The i-dependence of A/ c becomes stronger as n* is increased. The 
^-dependence of A/ c is closely related to the probability of finding a path ("saddle 
point" ) that connects the starting minimum to a different one. If such paths were 
relatively easy to find, then the transition probability would be insensitive to the 
value of t as long as it is not very short. If, however, paths to other minima are 
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few, a large number of configurations have to be explored before one of them is 
found. The t-dependence of A/ c would then be more pronounced and extend to 
larger values of t. To make the idea more concrete, we ignore the short-range time 
correlations among the configurations generated in a MC run and assume that they 
represent t independent samplings of configurations with free energy less than F Q +AF. 
Neglecting the rare return to the basin of attraction of the starting minimum after 
a transition to a different basin of attraction, the transition probability may then be 
estimated as P(n*,Af,t) = 1 - [1 -p(n*, A/)]* ~ l-exp(-tp), where p{n* , A/) < 1 
is the probability that a randomly chosen configuration with (3F < (3Fq + NAf 
belongs in the basin of attraction of a different minimum. One expects p to be zero 
if A/ < A/o(n*) where NUbTAJq is the height of the lowest free energy barrier, and 
p = g(n*,Af — A/o) for A/ > A/o where g(n*,x) grows continuously from zero as 
x is increased from zero. Combining this with the definition of A/ c , we obtain the 
relation g(n*, Af c (n*, t) - A/ (n*)) = ln2/t. Since A/ C (n*, h) - A/ C (n*, i 2 ) for fixed 
t\ < t2 increases with n* , the function g(n*,x) decreases (i.e. the difficulty of finding 
paths to other minima increases) as n* is increased at fixed x. 
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Figure 3. (a) Plots of A/ c , obtained for a L = 15 minimum, against 
(t/1000) -0 ' 35 for four values of n* . The dashed lines are the best straight-line 
fits, (b) The dependence of A/ C (n*,t = 15000) on n* for L = 15. The dashed 
line shows the best fit to a Vogel-Fulcher form (see text). 

The observed t-dependence of A/ c for all values of n* and all the minima in our 
study is well-represented by 

A/ c (nV) = Af (n*) + c(n*)t~ a , (3) 

with a in the range 0.25 — 0.40. Fits to this form with a — 0.35 for a minimum with 
L = 15 are shown in Fig. ||a. The values of A/o obtained from such fits with a fixed 
value of a are nearly independent of n*, but exhibit a dependence on the value of 
a, varying between and 0.5 for the L — 15 minimum. Similar results arc obtained 
for L = 12, with values of A/o between 1.3 and 1.5. The quantity c(n*) increases 
with n*. These results correspond to g(n*,x) ~ A(n*)x 1 ^ a with A(n*) decreasing 
with increasing n*. We conclude that the growth of the effective barrier height with 
increasing n* is primarily due to an entropic mechanism associated with an increase 
of the difficulty in finding low-lying saddle points that connect different glassy local 
minimum of the free energy. 
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The dependence of A/ c on n* is consistent with the Vogel-Fulcher law [[l8] which 
assumes the following form Jl9| for our system: 

Af c (n*) = a + b/(n*-n*), (4) 

where a, b and n* are constants. The value of n* obtained from fits of our data for 
A/ C (n*,i) to Eq.(^) with fixed a is nearly independent of t. This is consistent with 
the form of Eq.(Q) if a = A/o , b oc t~ a , and c oc l/(n* — n*). A/o is indeed nearly 
independent of n*, and we find that the ^-dependence of b and the n*-dependence of 
c are in agreement with the other two conditions. For the L = 15 case, we can fit the 
data for A/ c at t = 15,000 to the form of Eq.(f§) with a = (A/ = 0). The best 
fit, shown in Fig. ||b, corresponds to n* = 1.225, very close to the expected random 
close packing density n* cp ~ 1.23. The best fit to the L — 12 data with a ~ 1.0 also 
yields a similar value of n*. We conclude that the observed growth of the effective 
barrier height is consistent with the Vogel-Fulcher form. The increase in the effective 
barrier height as n* is increased from 0.96 to 1.06 is about 25fcsT, corresponding to 
a growth of the characteristic time scale of about ten orders of magnitude. Thus, the 
range of time scales covered in our study is comparable to that used in Vogel-Fulcher 
fits of experimental data, and much wider than what can be achieved in standard MC 
or molecular dynamics simulations. 



4. Discussion 



We close with a discussion of the connections between our results and those of spin- 
glass- like theories || 0, ^5| of the structural glass transition (see Ref. 22 for a 



discussion of the relation of our work with other recent studies of the behavior of simple 
liquids near the glass transition). These theories are based on the similarity between 
the phenomenology of the structural glass transition in so- called "fragile" § liquids 
and the behavior found in a class of generalized mean-field spin glass models || [26| 
with infinite-range interactions. At high temperatures, the free energy of these 
mean-field models, expressed as a function of the single-site magnetizations, exhibits 
only one, "paramagnetic" , minimum. As T is lowered, an exponentially large 
number of non-trivial local minima come into existence at a temperature T^, where 
a "dynamic transition" , characterized by a breaking of ergodicity, occurs. This 
"dynamic transition" does not have any signature in the equilibrium behavior of the 
system. A thermodynamic phase transition occurs at a lower temperature T c . In 
the suggested analogy between these models and the structural glass transition, the 
paramagnetic minimum of the free energy is identified with that corresponding to 
the uniform liquid, and the role of the non-trivial local minima of the free energy 
is played by the glassy local minima. The analogue of the "dynamic transition" at 
Td is thought to be smeared out in liquids. It has been suggested ||, 0, |||] that Td 
should be identified with the "ideal glass transition" temperature of mode-coupling 
theories [ p7[ . The temperature T c is interpreted as the "Kauzmann temperature" p8| 
at which the difference in entropy between the supercooled liquid and the crystalline 
solid extrapolates to zero. The relaxation time of the supercooled liquid is supposed 
to diverge at this temperature. Heuristic arguments suggest that this divergence is of 
the Vogel-Fulcher form |Q ^5| . 

Our results qualitatively support this scenario. We find a characteristic density 
at which a large number of glassy minima of the free energy appear. We do not know 
whether the number of glassy minima depends exponentially on the sample volume. 
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The configurational entropy associated with these minima decreases with increasing 
density because the width of the band over which the free energy of these minima is 
distributed increases with density. We have also found evidence for a Vogel-Fulcher- 
type growth of relaxation times driven by an entropic mechanism. 

There are, however, certain differences between our findings and the predictions 
of spin-glass-like theories. In our earlier work Jl4|, we found that the free energy 
of a typical glassy minimum becomes lower than that of the uniform liquid as the 
density is increased slightly above the value (n* ~ 0.8) at which the minimum comes 
into existence. In particular, the free energies of the glassy minima are substantially 
lower than that of the uniform liquid one for n* near n* ~ 0.95, the crossover 
density JlJ, |l5| above which the dynamics is governed by transitions among glassy 
minima. This is different from the behavior found in the spin glass models. Our 
results for the distribution of the overlap between different minima are also somewhat 
different from those for the spin glass models. Some of these differences may be due 
to finite-size effects which are probably significant for the small samples considered 
here. Also, fluctuation effects, which are unimportant in mean-field models, may play 
an important role in our system. A careful investigation of these issues would be 
interesting. 
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